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Abstract. 

We study resonating- valence-bond (RVB) states on the square lattice of spins and of 
dimers, as well as SU (A^)- invariant states that interpolate between the two. These states 
are ground states of gapless models, although the 5^/(2) -invariant spin RVB state is also 
believed to be a gapped liquid in its spinful sector. We show that the gapless behavior 
in spin and dimer RVB states is qualitatively similar by studying the Renyi entropy for 
splitting a torus into two cylinders. We compute this exactly for dimers, showing it behaves 
similarly to the familiar one-dimensional log term, although not identically. We extend the 
exact computation to an effective theory believed to interpolate among these states. By 
numerical calculations for the SU{2) RVB state and its S'C/( A") -invariant generalizations, 
we provide further support for this belief. We also show how the entanglement entropy 
behaves qualitatively differently for different values of the Renyi index n, with large values 
of n proving a more sensitive probe here, by virtue of exhibiting a striking even/odd effect. 
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1. Introduction 

A vast amount of recent theoretical and experimental effort has been devoted to the study 
of and search for spin liquids. A spin-liquid phase typically exhibits exotic behavior such 
as quasiparticles with fractionalized quantum numbers, or spin-charge separation. Such 
behavior is manifestly non-perturbative, and so requires advanced theoretical tools to 
analyze. 

A vast amount of recent theoretical effort has been devoted to understanding how 
the entanglement properties of a particular state illuminate its physical properties. Since 
a major lesson of condensed matter physics in recent decades was that a phase often 
can be characterized by a model state, a feature of studying entanglement is that states 
can be studied directly, without needing to know a specific Hamiltonian. Since studying 
entanglement properties of a model state is often both analytically and numerically more 
tractable than analyzing a model Hamiltonian, much recent progress has been made in 
condensed matter physics from the former. 

The behavior of quantities such as the Renyi entropy thus can be used to explore spin- 
liquid states directly. The quantum spin liquids best characterized theoretically are gapped, 
where a gauge symmetry typically emerges. An important and useful concept applicable here 
is topological order, which relates the emergent gauge symmetry to a topological degeneracy, 
and the related concept of topological entropy [1, 2, 3]. However, many experimental systems 
currently identified as good candidates for spin liquids appear to be gapless. Compared to 
their gapped counterparts, these phases are less amenable to simple unifying attributes 
such as topological order. Although it is likely that other quantities, related perhaps to 
correlations or manifestations of long-range entanglement, can be used to characterize gapless 
phases, few viable examples of models - particularly those tractable by large-scale numerical 
simulation - are known to contain gapless spin liquids. 

We study a particular family of states, resonating-valence-bond (RVB) states on the 
square lattice, that are gapless and exhibit spin-liquid-type behavior [4, 5, 6]. Our main 
(but not only) tool is to analyze a subleading term in the Renyi entropy particularly useful 
for characterizing gapless systems [7] . The Renyi entanglement entropies Sn are defined as 

Sn = ^lnTvp^^, (1) 
I — n 

where pa is the reduced density matrix for a region A, and the von Neumann entropy Si is 
defined as the limit n ^ 1. The term we analyze describes the entanglement between two 
regions formed by cutting a torus into two cylinders. This geometry is useful for studying 
gapless systems because the length of the boundaries separating the cylinders is independent 
of the area of the cylinders. Because of the long-range correlations, the entanglement entropy 
of gapless systems will be sensitive to varying the size of the cylinders. This dependence was 
studied for three completely different gapless systems in [7], and argued to have universal 
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behavior. One of the results of this paper is to provide additional evidence in support of this 
assertion. 

One model with a gapless RVB ground state, the square-lattice quantum dimer model, 
is quite amenable to an analytic treatment. For example, all equal-time correlators diagonal 
in the dimer basis are those of a classical dimer model [4]. These can be computed exactly 
directly in the dimer model by using Pfafhan techniques [8, 9, 10], or in the continuum 
limit by utilizing a two-dimensional classical free-boson field theory [11, 12]. Such exact 
computations have been done for entanglement quantities as well [13, 14, 15, 16, 17, 18]. In 
this paper we extend these methods to compute the two-cylinder entanglement entropy for 
the quantum dimer model exactly in the scaling limit. For Renyi index n greater than a 
specific (non-universal) value nc, it is simply given by a ratio of cylindrical partition functions. 
From this we find for example that the striking even-odd effect observed numerically in [7] for 
the S'[/(2)-invariant RVB state was not a finite-size artifact; it also occurs for the quantum 
dimer model, and our computation illuminates its origin. 

One lesson from our results is that universal entanglement properties can be qualitatively 
different for values of the Renyi index n. For n > ric^ the two-cylinder entanglement entropy 
depends strongly on whether the cylinders are of odd or even length. This is not shocking 
from the point of view of the quantum dimer state, because the correlators exhibit similar 
behavior. However, at the von Neumann point (n ^ 1), no such even-odd effect is possible. 
Indeed, the shape of the entanglement entropy has to be a concave-down function of as 
required by the strong subadditivity property [19]. Other values of n are not constrained 
by subadditivity, and we will derive this even/odd effect explicitly for the dimer RVB state. 
Thus only for n > ric is the Renyi entropy sensitive enough to include this effect. This 
effect, whereby a Renyi entropy may exhibit qualitatively different behavior for some higher 
n, is reminiscent of the manifestation of finite-temperature criticality in Renyi entropies for 
n > 1, but not n = 1 [20]. 

It is possible to interpolate between the 5'[/(2)-invariant RVB state and the quantum 
dimer state by considering 5'[/( A^)-invariant generalizations of the former; many aspects 
of the quantum dimer state then can be extracted from the ^ oc limit. Such a 
correspondence can be pushed further: Ref. [21] defined a cluster expansion of the RVB 
"loop gas" [22] in terms of interacting dimers in the large limit. One recovers to first non- 
trivial order the interacting generalization of the quantum dimer model studied in [23, 24]. 
Most strikingly, the critical exponent a extrapolated to = 2 turned out to be in good 
agreement with the previous numerical studies [5, 6]. This expansion also provides some 
insight in the study of analogous RVB states in three-dimension [25, 26]. 

Whereas one cannot do exact computations directly in these models outside of the dimer 
case, the field theory computation can be extended to a full line of critical points by varying 
the stiffness of the classical boson. This line of critical points likely describes the scaling limit 
of the interacting quantum dimer model [23, 24, 21]. We expect that the SU{N) version 
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studied here also will have the same universal behavior, and by comparing our quantum 
Monte Carlo results to the expansion developed in [21], we provide strong evidence that the 
scaling limit of the = 2, 3, 4, 5 theories indeed are points on this line. Moreover, we provide 
consistency checks that the two-cylinder entanglement entropy also behaves universally for 
= 2, 3, 4. Our results therefore strongly imply that at least the gapless part of the square- 
lattice RVB state behaves qualitatively similar to the quantum dimer model; there is no 
phase transition as is varied. 

The paper is organized as follows: In section 2, we review the definition and important 
properties of the RVB states. We derive our general result for the Renyi entanglement 
entropy of many critical Rokhsar-Kivelson type states, including the quantum dimer state, 
in section 3. Namely, for n > ric the Renyi entropy is given by a simple ratio of 2d classical 
partition functions. Section 4 explores the consequences of this result in the simplest case of 
the square lattice quantum dimer state. We also discuss the emergence of a strong even-odd 
effect, similar to that first observed in [7]. We then turn our attention in section 5 to the 
SU{N) RVB state. By quantum Monte Carlo simulations, we show that while the spin-spin 
correlations decay exponentially, the dimer-dimer correlations are critical. The exponent 
approaches the dimer exponent {a^ = 2) SlS N increases, and agrees well with a previous 
large computation [21] and earlier numerical results at A^ = 2 [5, 6]. We also present 
further evidence of the underlying Coulomb-gas structure for all A^. Finally, in section 6 
we present numerical simulations of entanglement for the SU{N) RVB wave function. The 
behavior qualitatively resembles that of the dimer RVB state, and therefore provides strong 
evidence of universal behavior, in agreement with our main result of sec. 3. 

2. The RVB states and their equal-time correlators 

A paradigm for a spin liquid is the two-dimensional resonating-valence-bond (RVB) state 
[27]. Much effort has gone into its study, since the classic work in the late '80s [28, 22]. In 
this section we introduce and review some of the known properties of these states. 

2.1. The SU{2) RVB state 

A valence bond is a spin singlet, and a valence-bond state is one where each spin forms 
a valence bond (i.e. a singlet) with one other spin. The nearest-neighbor SU{2) RVB 
wavefunction is the equal-amplitude superposition of all nearest-neighbor valence-bond (or 
singlet) states of spin-1/2 particles fixed at the sites of some lattice. For the A/'-site square 
lattice, the SU{2) RVB state is 




(2) 
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so each spin i on one sublattice is in a singlet with one of its nearest neighbors ji on the other 
sublattice. A valence bond between neighbors i and j can be labeled by a dimer on this link 
between i and and so each valence-bond state can be viewed as a dimer configuration C 
with exactly one dimer touching each site. Our computations probe a state directly and so 
do not require a Hamiltonian. However, it is worth noting that this SU (2) RVB state is the 
exact ground state of a quantum spin Hamiltonian with local interactions [29] . 

This RVB state breaks no symmetry, and so is a natural candidate for a state exhibiting 
spin-liquid behavior. For example, if the two-dimensional square lattice is periodic (i.e. is 
a torus) and has an even number of sites in both directions, there is a winding number for 
each direction. These are defined simply by counting the dimers around a closed path going 
straight around a cycle, weighting each dimer by +1 for an even number of steps along the 
path, —1 for an odd number of steps, with no contribution from empty links. An RVB state 
can be defined for each value of the winding numbers, and each will be the ground state 
of a local Hamiltonian. The number of ground states will then depend on the number of 
non-trivial cycles, a basic characteristic of topological order. 

For the square lattice, extensive numerical analysis indicates that dimer-dimer 
correlators in the RVB state decay algebraically [5,6]. A theorem of Hastings guarantees that 
when any local operators have an algebraically decaying correlator, any local Hamiltonian 
with this as a ground state must be gapless [30] . Thus the square-lattice RVB state could only 
describe a gapless system. However, it has long been known, both from numerics and strong 
analytic arguments, that spin-spin correlators are exponentially decaying [28]. Coupled with 
the fact that no local order parameter has been found [5, 6], it seems very possible that there 
is a spin gap, and that the state effectively behaves as a spin liquid. In this paper we do not 
address this question directly. Instead we focus on the gapless behavior of nearest-neighbor 
RVB states on the square lattice. 

2.2. The dimer RVB state 

One key tool is to study the dimer RVB state [4]. Here one essentially forgets the underlying 
spins: the only degrees of freedom are the nearest-neighbor dimers. Each configuration of 
dimers is a linearly independent basis element of the Hilbert space of the quantum dimer 
model, and these dimer "states" span the Hilbert space. As above, a dimer configuration C 
has exactly one dimer touching each site of the lattice, i.e. the dimers are close packed and 
hard-core. The RVB state for dimers is the equal-amplitude sum over all these states: 




(4) 
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Different dimer configurations are defined to be orthogonal: 
{Dc\Dc') = 5c ' 



(5) 



The "RK" Hamiltonian [4] is a simple local Hamiltonian with this dimer RVB state as its 
ground state. 

One very convenient feature of the dimer RVB state is that any correlator diagonal in the 
dimer basis can be computed exactly, because they are exactly those of the two-dimensional 
classical dimer model. This follows immediately from the definitions of the state and the 
inner product. For example, the normalization of the dimer RVB state is the partition 
function of the classical dimer model, the number of distinct dimer configurations: 



which behaves asymptotically on the square lattice as ^ (1.79162 . . Y [8]- A classic result 
of classical statistical mechanics is that any correlator in the classical dimer model on any 
planar lattice can be expressed as a Pfaffian [8, 9]. This is equivalent to saying that the model 
can be rewritten in terms of free fermions [31]. For the dimer-dimer correlation function, 
the asymptotic behavior can easily be obtained using this expression. For the square lattice 
long ago this was shown to be algebraically decaying: 



with exponent a = 2 [10]. The coefficient C depends on both the direction and the parity 
of the number of steps along any path connecting the two dimers; in some special cases 
it vanishes and so subleading terms dominate. In fermionic language, the dimer creation 
operator is a fermion bilinear, and so is of dimension 1. Algebraic decay occurs only for 
bipartite lattices; for non-bipartite lattices, correlators are exponentially decaying, a fact 
crucial in showing that the triangular-lattice quantum dimer model is in a gapped RVB 
phase, i.e. has topological order [32]. 

It is also possible to generalize this construction, with "enforced" orthogonality (5), to 
other models. For example starting from a classical six-vertex model, we get a quantum six- 
vertex wave function [33]. This wave function also has algebraically decaying correlations, 
with a continuously varying exponent a as a function of the vertex weights on its critical 
line. The classical model, although in general more complicated than the dimer model, is 
nevertheless integrable. It is also constructed with a similar set of constraints, the fully 
packed and hardcore constraints being replaced by the ice rule. All the analytical results we 
will obtain for the entropy of quantum dimers can be straightforwardly generalized to the 
quantum six-vertex case. 

The inner product is an important difference between the dimer and SU (2) RVB states. 
In the RVB state (2,3), the inner product is the usual one for spins: each different spin state 
in the diagonal basis is orthogonal to the others. As is easily shown [22], this means that 
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different valence-bond states in the spin model are not orthogonal. There is a convenient 
loop-gas description to represent their inner product. One places the two dimer graphs on 
top of each other, obtaining the "transition" graph. Because of the requirement that there 
be exactly one dimer per site, this means that the transition graph is comprised of closed 
loops. The shortest loops are of length two, when both graphs have a dimer on the same 
link. Using the spin inner product then gives [34] , 

{Vc\Vc>) = 2^^-^/\ (8) 

Til being the number of loops in the transition graph. Despite this overlap between dimer 
configurations making analytic computations much more difficult, its presence is convenient 
in the algorithm used for computing the second Renyi entropy using quantum Monte Carlo 
simulations [35], since "swapped" basis configurations (described in Section 6.1) will always 
contribute a non-zero expectation value, even with a naive QMC updating algorithm. 

Despite this important difference between the dimer and SU{2) RVB states, their 
obvious similarity makes it plausible that the dimer-dimer correlator will behave similarly in 
the two cases. Further evidence for this fact is that rewritten in loop language, the spin-spin 
correlator corresponds to the probability that the two spins are on the same loop. This 
decays exponentially for loop models with weight per loop greater than two [36], and the 
inner product (8) amounts to a weight per loop of 4 [28]. This exponential decay arises 
because the classical partition function is dominated by configurations with many short 
loops. Thus correlators in the SU{2) RVB state are dominated by short-loop configurations, 
i.e. the dimers. Whereas the presence of longer loops renormalizes the exponent a, it seems 
very likely that dimer-dimer correlators in the RVB state will decay algebraically like they 
do for the dimer RVB state. 

This expectation was convincingly demonstrated numerically [5, 6]. The dimer-dimer 
correlator for dimers was found to be of the form (7) with exponent close to a 1.2. This is 
an even weaker falloff than in the dimer RVB state, so the presence of the loops in the inner 
product of the SU (2) RVB state only makes this correlator longer-ranged. By the theorem 
of [30] , algebraic decay of local operators in a two-dimensional ground state requires that any 
local Hamiltonian with this ground state be gapless. Thus not only is the quantum dimer 
model with RK Hamiltonian gapless, the spin Hamiltonian of Ref. [29], whose exact ground 
state is the SU (2) RVB state, must be as well. 

2.3. The SU{N) RVB state 

It is possible to interpolate between the dimer and SU{2) RVB states, either by dressing 
the lattice itself [37], or by modifying the SU{2) symmetry to SU{N). In the SU{N) RVB 
state, the "spins" have = 25^ + 1 components, and an SU (N) singlet between sites i and 
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j is given by 

[^,j1 = ^^i=f E (-ir-"H,®|-m), . (9) 

When i and j are nearest neighbors, such a singlet can be viewed as a dimer occupying the 
hnk between i and as in the SU{2) case. The nearest-neighbor RVB wave function is then 
defined as an equaUy weighted superposition of these valence-bond states (2) as before, with 

\^c) = ^,\{[^.3^] . (10) 

i=l 

For SU{N) valence-bond states, the orthogonality relation becomes 

{Vc\Va) = ±{2S + ir-^'\ (11) 

where ni is again the number of loops in the transition graph formed by superimposing the 
two dimer configurations C and C . Notice that ni can be at most A/"/ 2, and this happens 
only when all the loops are trivially fiat, i.e when the two dimer configurations are identical. 
We have however 

{Vc\Va) = 5c,c' (12) 

in the limit ^ oc (or equivalently oc). As a result, any correlation calculated in 

the limit N ^ oo will be identical to that in the corresponding dimer model. It is however 
important to stress that this argument strictly only holds at the level of correlation functions. 

We will present in this paper a variety of results for the SU{N) case, and show that 
there is no sign of a phase transition as is varied, thus giving further evidence for the 
similarity in behavior between the SU (2) and dimer RVB states. We will see, for example, 
that for A" = 2,3,4,5 the dimer-dimer correlator is algebraically decaying, with exponent 
depending on A^. 

Given the changing exponents, the SU{N) RVB states are obviously not identical to 
that for dimers, even in the scaling limit. In fact, we will present strong evidence that equal- 
time correlators for the dimer and SU (N) RVB states can all be written in terms of a free 
boson, or Coulomb gas. 



2.4' The Coulomb- gas approach to classical dimers 

A well established and very useful fact about the classical dimer model is that long-distance 
properties of the dimers are described by a free bosonic field theory, or Coulomb gas. This 
can be seen either by taking the continuum limit of the fermionic description (see e.g. [38]) 
and then bosonizing the fermions, or by rewriting the dimers in terms of a discrete "height" 
degree of freedom [39, 40, 12]. We follow the latter approach, because a careful treatment of 
the boundary conditions is necessary to obtain our results for the two-cylinder Renyi entropy. 
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Figure 1. Illustration of the height shift between the upper and lower boundaries in the 
odd case. Left: Even Ly case. The coarse-grained height difference between top and bottom 
is Ah = (htop) — (/^'bottom) = ('7r/2)(9/2 — 1/2) = 27r. Right: Odd Ly case. We have 
Ah = (7r/2)(13/2 - 1/2) = Stt. 



The first step in the mapping of classical dimers on to the Coulomb gas is to associate 
an integer with each face of the lattice, as illustrated in fig. 1. These integers are defined 
by the following construction. Set value on one face (say the lower leftmost one) to be zero. 
Then, turning counterclockwise around a site of the even (resp. odd) sublattice represented 
by black squares (resp. white squares) the integer changes by —3 (resp. +3) if it crosses a 
dimer, +1 (resp. —1) otherwise. Since by construction these integers are different on adjacent 
sites, in the bulk one defines a smoother height variable h{i^j) on the sites of the lattice, 
by averaging the values of the integers on four faces surrounding each point (and then for 
convenience's sake multiplying by 7r/2). Then the states where the dimers are arranged in 
columns correspond to a constant value of the height. 

An important subtlety in this mapping is the boundary conditions, which will play a 
crucial role in understanding the even/odd effect in the two-cylinder Renyi entropy. With 
periodic boundary conditions and an even number of sites around this cycle (as in the left 
and right sides of the figure 1), note that means that heights can jump by 27r times an integer 
on from one side to the other. With open boundary conditions, as we have at the ends of the 
cylinder (the top and bottom of the figure), the natural definition of the boundary height 
comes from averaging the integers along the boundary. It is then easy to check that for 
a cylinder of even circumference and length Ly^ the height difference between top and 
bottom satisfies 

Even/odd effects have long been studied in the dimer model [41, 42, 43], but to our knowledge 
they have never been interpreted in the context of the height mapping. 

The standard Coulomb-gas arguments [36] then imply upon coarse-graining the height 



10 



variable h{i^j) into a continuous field, the scaling limit of the classical dimer model is 
described by a free scalar field h{x^y) with Euclidean action 



The coupling is often known as the "stiflFness"; we show below that for the dimer case 
it is set to /^D = 1/2. This free-boson field theory is one of the simplest examples of a 
two-dimensional conformal field theory [44]. It is ubiquitous in two-dimensional statistical 
models, as well as in condensed matter, where it describes Luttinger liquids. In order to be 
consistent with periodic boundary conditions, field values shifted by 27r must be identified; 
in conventional language it is said that the field h is compactified on a circle of radius 1: 
hr^ h + 271. 

The dimers are viewed as elementary electric charges, i.e. are created by the field cos(/i) 
(the simplest function of h that satisfies the compactification h ^ h + 27r). The standard 
computation [36] gives for the dimer-dimer correlator decay exponent: 

a = -. (15) 

This yields the stiffness for the dimer model k^d = 1/2- Elementary magnetic charges can 
also identified with monomers, and the corresponding exponent describing their two-point 
function is given by /3 = so /3d = 1/2. Once the stiffness fixed to the appropriate value, 
Eq. (14) describes all the universal properties of the classical dimer model. By virtue of the 
correspondence between correlators in the classical model and those in the dimer RVB state, 
the same exponents apply to the correlators in the latter as well. The same Coulomb-gas 
techniques can be used on the six-vertex model as well, yielding a value of continuously 
varying with the couplings [36] . 

In this framework the quantum dimer and quantum six- vertex models can also be written 
as in terms of a free scalar field in the 2+1 dimensional "quantum Lifshitz" model [12, 33]. 
This is a field-theory version of an RK Hamiltonian, where the equal-time correlators of the 
2d quantum model are those of a 2d conformal field theory. This model is critical but not 
Lorentz-invariant, with dynamical critical exponent z = 2. 

Because of the arguments above, it is reasonable to hope that all the equal-time 
correlators in the SU (N) RVB state for any can be described by this Coulomb gas, albeit 
with stiffness depending on A^. Substantial evidence for this was provided in the SU (2) case 
in [6], where several universal quantities, including the dimer-dimer and monomer-monomer 
decay exponents, were found numerically. With appropriate identification of the operators, 
each yields a independent value for the stiffness. The numerics gave approximately the same 
value of the stiffness k.2 ^ -83 for each, providing strong evidence that the Coulomb gas 
description applies to correlators of spin singlets in the SU (2) RVB wave function. 




(14) 
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3. Two-cylinder Renyi entropy from partition functions 



In this section, we discuss the shape-dependent subleading contribution:]: to the entanglement 
entropy of two-dimensional gapless systems discussed in [7] . For Renyi index n larger than a 
certain critical value nc, we derive it exactly for the scaling limit of the square-lattice dimer 
RVB state, the quantum six- vertex state, and the entire quantum critical line of the quantum 
Lifshitz theory Thus we expect it to apply to the scaling limit of the SU (N) RVB state as 
well, with an appropriate parameterization of in terms of A^. 

In particular, we show that when a torus is split into two cylinders labeled A and 5, 
the Renyi entropy describing the entanglement between the two cylinders can be rewritten 
in terms of conformal field theory partition functions Z as: 



Below we detail the precise definition of the CFT partition functions, but the important fact 
is that for the cases at hand, these are known. Thus this expression gives an explicit formula 
for this entropy, exact in the scaling limit. The value of ric depends on the specifics of the 
model, but ric < 2 in all the examples we will consider, so that Eq. (16) applies to the second 
Renyi entropy to be studied numerically in later sections. Although our derivation applies 
only to the quantum Lifshitz ground state (i.e. RVB-type states that can be written in terms 
of a two-dimensional classical free boson), we believe it possible that an analogous formula 
will apply more generally, in particular to any theory with an RK-type Hamiltonian. 

3.1. Entanglement entropy as a Shannon entropy 

We will now detail how to derive Eq. (16) in the simplest case of quantum dimers on the 
square lattice. This requires a slight generalization of two ingredients: the mapping to a 
classical Shannon entropy of Ref. [15] (see also [45]) and the boundary phase transition 
argument of Ref. [18] (see Sec. 3.2). Although this derivation does not apply to the SU{N) 
RVB case, we will later explain why we think Eq. (16) should still apply. 

The orthogonality of the different dimer configurations in the dimer RVB state allows 
for huge technical simplifications. As was shown in Ref. [15], the entanglement entropy can 
be expressed as a classical Shannon entropy. We give here only a brief summary of the 
results, and refer to [15] for more details. 

Let us cut our torus in two subsystems A and 5, as is shown for example in Fig. 2. 
We label the dimer configurations along the boundaries as \a) and The main difference 
with [15] is that there are two boundaries here, but the the same arguments apply. The von 

i Of course, the non-universal leading "area law" contribution is shape-independent. Possible additional 
universal logarithmic terms, not encountered in the two-cylinder geometry, may also be shape-independent. 
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Figure 2. Bipartition of the 8x8 torus. The two boundaries between A and B are 
emphasized by thick dashed hues. Links crossing the lower boundary are defined to be in 
A, whereas those crossing the upper boundary in B. 



Neumann entropy can be recast as a classical Shannon entropy 

S = -^Pa,f,^^Pa,f, , (17) 

where the Pa^,^ are the probabilities of a given boundary configuration |cr),|/x) between A and 
B. The Pcr^i^ are precisely the eigenvalues of the reduced density matrix. Said differently, the 
Schmidt decomposition of the ground state reads 

i*i^) = Ev^iOiO' (18) 



'A 



where IV^^^) is a superposition of all dimer states in A compatible with the boundary 
conditions a and /x, and the same goes for B. All these vectors are mutually orthogonal: 

{i^^,M'y)=S.a'S,,>Snn' • (19) 
This mapping therefore allows the Renyi entanglement entropy to be rewritten as 

Sn = Y^l^ [J2^.,A . (20) 

We name this quantity the Renyi-Shannon entropy in the following. The probabilities can 
be conveniently rewritten in terms of dimer partition functions as 

Pa,, = (21) 

where Z^r,^ is the number of dimer configurations compatible with the boundary configuration 
I a), l/x), while Z is the number of dimer coverings on the whole torus. 

It is important to emphasize the main reasons why such a remarkable quantum-classical 
correspondence holds. While (18) remains true for all the states we consider in this paper, 
the orthogonality (19) is guaranteed only if three other conditions are satisfied: 
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(i) The orthogonality of quantum states corresponding to different classical configurations. 

(ii) A certain type of local constraint [15], which for Z2 variables on the bonds of the square 
lattice (e.g. dimer occupancies) amounts to: if three of the variables around each site 
are specified, then the value of the fourth is uniquely determined. 

(iii) Interactions are only between bond variables sharing a common site. 

Quantum dimers obviously satisfy both these requirements. The same goes for the quantum 
six- vertex state, (ii) being guaranteed by the ice rule. The SU{N) RVB states satisfy (ii) 
and (iii) but not (i) (see (8)). An example of a Coulomb-gas wave function that satisfies 
(i) and (ii) but not (iii) is a quantum dimer wave function built out of an interacting dimer 
model [23, 24]. 

Returning to the dimer RVB state, these expressions make it possible to use the 
computer to find essentially exact values for the Renyi entropy. Using free-fermion techniques 
described in Appendix A, each probability can be expressed as a product of determinants 
of two matrices of size ^ L/2, and computed numerically in a time of order ^ L^. Using 
translational symmetry along the x axis as well as the conservation of winding numbers, 
we are able to compute any Renyi entropy up to machine precision in a time of order 
^3/2 ^ rpj^^ results for a 20 x 20 torus are summarized in Fig. 3, where the two cylinders 
are of length £y and Ly — £y. The figure makes it obvious that the Renyi entropy exhibits 




6 I ^ ^ ^ J 

5 10 15 20 

£y 

Figure 3. Two-cylinder Renyi entanglement entropies 5'i/2, Si^ S2 and for the dimer 
RVB state on the square lattice in the torus geometry. Here = Ly = 20. 

completely diflFerent behavior for diflFerent values of n. In the "replica" phase n < 1, the 
curves are relatively fiat. In the thermodynamic limit we expect them to go to universal 
functions, which in principle may be computed within boundary CFT, although we will not 
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study this phase here. For n large enough, we observe an even-odd effect hke that observed 
for the SU{2) RVB state in [7]. We explain the origin of this transition in dimer RVB state 
in the next subsection. 



3.2. The entropy in the locked phase n > Uc 

Here we derive the relation of the two-cylinder entanglement entropy of the dimer RVB state 
to CFT partition functions in (16) for n > Uc- This follows fairly simply from the description 
of the long-distance behavior of classical dimers in terms of the free-boson action (14). 

When one utilizes an effective field-theory action to describe the scaling limit of a lattice 
model, one must include all terms in the action consistent with the symmetries of the field 
theory. As is familiar from the derivation of the Kosterlitz-Thouless transition in field theory, 
operators of the form 

Vd = - cos {d{h- ho)) (22) 

with d integer are consistent with the compactification of the boson h ^ h + 27t. If such 
a term is relevant, then the action flows to a configuration where the height h is locked to 
one of the minimal values h = ho mod 27r/rf. For dimers on the square lattice, it has long 
been known that any such term allowed are irrelevant in the bulk; this is why the action (14) 
applies and the correlators are algebraically decaying. However, we showed in the preceding 
subsection 3.1 that the probabilities used to compute the Renyi entropy depend on the 
boundary conditions. Thus the crucial observation of [18] (see also [15]) is that one needs to 
include terms of the form (22) localized at the boundary. 

To be precise, in the continuum limit the probabilities p^^i^ depend now on the 
configurations of the field p((/>). The probability of observing field configuration (j) at the 
boundary is still Gaussian 

p{<P)^exp{-SM), (23) 
and raising this probability to the n— th power yields 

[p{4>)r = exp{-nSM) = exp(-^„,[0]). (24) 

Therefore, the system near the boundary effectively feels a stiffness k.^ = nn. Standard 
Coulomb-gas/field-theory computations (see e.g. [44, 36]) then give the boundary scaling 
dimension of (22) to be d^ /{2nn). As a consequence, Vd is irrelevant as long as d^ > 2nK. 
Thinking of n as a continuous parameter, this defines two regions, separated by a (boundary) 
Kosterlitz-Thouless phase transition. The critical value of n in the computation of Sn is 
therefore 
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where rfmin is the smaUest d allowed by the lattice symmetries. For dimers on the square 
lattice rfmin = 1 and = 1/2, so that the critical value is 

ric = 1. (26) 

The two qualitatively different types of behavior seen in fig. 3 indeed occur on opposite 
sides of Tie- It is important to note that even though Sn should take an universal form in 
both phases, the value itself of the critical Renyi parameter depends on a degeneracy rfmin, 
which is non-universal. For example, rfmin = 2 for the quantum six-vertex wave function, 
and dmin = 3 for dimers on the hexagonal lattice. 

This observation allows us to derive the expression for Sn valid for n > ric- In this 
phase, the boundary operator is relevant. This means that here the boundary value of the 
field locks on to a minimum of Vi. This minimum must be fixed along the entire boundary; 
this is commonly known as Dirichlet boundary conditions on the boson. We refer therefore 
to the phase with n > ric as the "locked" phase. 

A constant value of the field along the boundary means that the action is at its minimum 
there, so this is the maximum value of the probability Pmax- Universal contributions to the 
Renyi entropy coming from the sum (20) are therefore dominated by the configuration: 

Tl 

Sn>nc ^ In(pmax) (27) 

I — n 

This result implies that in the scaling limit there is an "entanglement gap" between Pmax 
and the other values. We have confirmed directly the presence of the entanglement gap 
for the dimer RVB state by numerical evaluation of the probabilities using the free-fermion 
techniques detailed in Appendix A. This is distinct from the quantum Hall case [46], where 
there is no entanglement gap, in correspondence with the gapless edge excitations there. 

In terms of dimers, the configuration with maximum probability corresponds to no 
dimers crossing the boundaries between A and B (see Fig. 2). This is a huge technical 
simplification, since the probability of having no dimers across the cut factorizes into pieces 
coming from region A and from region B. Namely, define Zcy\{Lx^ Ly) to be the partition 
function (the number of dimer configurations) for a cylindrical region of size Lx in the periodic 
direction and Ly in the other, with boundary conditions corresponding to no dimers leaving 
the region. Likewise, Ztorus{Lx^ Ly) is the number of dimer configurations with periodic 
boundary conditions. Then it follows from (21) that when an Lx x Ly torus is split into 
cylinders of lengths £y and Ly — £y^ 

Zcy\{Lx^ iy)Zcy\{Lx^ Ly — £y) 

Pmax ^ TT • i^oj 

The analogous formula still holds if space is a cylinder instead of a a torus, simply by 
replacing the denominator with the appropriate cylinder partition function. Note that the 
non-universal bulk parts of the partition functions cancel in this expression. The non- 
universal boundary parts do not, and contribute to the area law in (27). Thus putting 
(28) together with (27) and filtering out the area law term gives (16). 
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For the dimer RVB state, the partition functions on the cyhnder can be computed 
exactly both on the lattice and in the continuum. The lattice dimer result is detailed in 
Appendix C. For the free-boson field theory, these partition functions are known explicitly 
[47]; we will discuss their application to the two-cylinder Renyi entropy in the next section 4. 
In particular, we will show how the even-odd effect is a signature of this locked phase n > ric, 
arising from the different boundary conditions necessary for the even and odd sectors. 

It should also be possible to extend these results to the region n < ric, using boundary 
CFT methods [15, 16, 17, 18]. Here, however, the boundary conditions will not allow the 
factorization of the probability into pieces coming from the two regions. Oshikawa[16] has 
already treated the y = £y/Ly = 1/2 case for the closely related cylinder geometry, using the 
replica approach n G N. There will possibly be subtleties involving analytic continuation 
in n, similar to the two-interval Id calculation [48, 49]. The study of the "marginal" case 
n = Tie would presumably be even more challenging [50, 18]. 

4. Entanglement in the dimer RVB state 

This section is devoted to the exploration of the consequences of Eq. (16) for the two-cylinder 
Renyi entropy of the dimer RVB state. We will calculate explicitly the partition functions, 
and so find results exact in the scaling limit. 

The main subtlety in this computation is how to account in the field theory for the 
distinct results occurring when the cylinders are of even and of odd length, as apparent 
in the numerical results in fig. 3 for the dimer RVB state and those in [7] for the SU{2) 
RVB state. The reason this is possible in the field theory is that a fixed/Dirichlet boundary 
condition is in fact a family of boundary conditions, corresponding to the specific value /iq 
that the field takes at a boundary. This value can be always be shifted overall, since the bulk 
action used to compute the partition functions is independent of it. However, on a cylinder, 
there are two separate boundaries, and what cannot be shifted away is the difference of the 
values on the two ends [51, 52]. The general arguments given above in section (3.2) do not 
specify this difference; it follows rather from the particular microscopic model. For dimers, 
we showed in section 2.4 that 

A/i = 2^{w + a) , ^ G Z, (29) 

where a = for even Ly and a = 1/2 for odd Ly. 

The computation of the free-boson partition function is completely standard; see e.g. 
[53, 54, 47]. Since the action is quadratic in the bosonic field, it can be decoupled into an 
oscillator part and a classical part. Only the classical part affected by the compactification 
/i ^ /i + 27r, and by the boundary conditions. The partition function for Dirichlet boundary 
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Figure 4. Numerical results for the second Renyi entropy S2 as a function of the subsystem- 
ratio y = iy/Ly for Lx = 20. Black dots are the data for square-lattice dimer RVB state, 
while the CFT predictions are given by Eq. (33) in the even case and Eq. (34) in the odd 
case. Left: Torus with aspect ratio Ly/Lx = 1. Right: Torus with aspect ratio Ly/Lx = 2. 



conditions on both ends of the cylinder, with the height differing by 27ramod27r, is given by 

CO +CO 

w=l w=—oo 

It is convenient to rewrite such sums in terms of the Jacobi theta and Dedekind eta functions 
defined in Appendix B, to take advantage of many elegant identities. The even case a = 
can be written, using (B.2,B.7,B.11,B.14,B.16), as 



^DDf^ ^3(2r) 

^cyi (r) = ^ , (31) 



where we adopt the conventional variable r = iLy/L^, and have set k = 1/2, the dimer 
value. For the odd case a = 1/2, we have, using this time (B.2,B.5,B.11,B.15,B.16): 



= ^ • ^^^^ 



Plugging these results into Eq. (16) and using the torus partition function [44, 47] gives for 
the two-cylinder entanglement entropies 

.(even) . _ ^ 1^ ( ^^^f ^ ^3(2yr)W^\ . .33. 

" ^ " 1 - n W^r)e,{r/2) r7(2yr)^(2(l - ^)r) ) ' ^^^^ 

" ^ " 1 - n W'^r)e,{Tl2) r7(2yr)r)(2(l - ^)r) ) ' ^^^^ 
Fig. 4 shows a numerical test of the universal shape Eq. (33) for the dimer RVB states with 
two different aspect ratios Ly/L^ = 1 and Ly/Lx = 2. Obviously, these n > ric results are 
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quite different in the even and odd sectors. The difference is proportional to n/{n — 1), and 
therefore shghtly diminishes with increasing n, but does not go away. 

This even/odd effect disappears as Ly/L^ is increased; it foUows from the CFT 
expressions that this occurs exponentiaUy quickly. It does not approach the celebrated 
Id result[55] oc lnsin(7r^) in the effective Id limit Ly/ oc, because the RK Hamiltonian 
becomes gapped with correlation length ^ (x m this limit. 




0.4 0.6 

y = (-y/Ly 



Figure 5. Comparison of numerics and analytics in the exactly solvable case 5*00 • Black 
points are the numerical data for all winding sectors, and black crosses represents the 
numerical data in the W = (0, 0) sector. Left: = 16, the data is compared with Eqs. (33), 
(34), (35) and (36) up to a single constant. Right: Same procedure for = 96. 



The relatively large finite-size effects apparent in the plots should disappear in the 
thermodynamic limit. To illustrate the slowness of the convergence, we study the exactly 
solvable case n ^ oc, where the agreement between CFT and the lattice can be made 
rigorous. It follows from (20) that only the largest probability (28) contributes: = 
— Inpmax, and so numerical analysis is much easier. We checked system sizes up to Lx = 704, 
and found that while indeed there are strong finite-size effects, the data converges in the end 
to the CFT result. The slow convergence is illustrated by the numerical and analytical results 
in performed in Fig. 5, for the aspect ratio Ly/L^ = 1 and the two system sizes = 16 
and Lx = 96. We also computed s^o in the W = (0,0) winding sector directly for the dimer 
model using the techniques described in the appendices, finding 

7^{Tyes{2y{l - y)T) \ 

7^{2yTU2{l - y)T) J ' ^ ' 

il{2Tfe,i2yil - y)T) \ 

7^{2yTM2{l - y)T) J' ^ ' 

Even though the winding numbers are zero, the result still depends on the compact ification 
radius because of winding fluctuations at the boundaries. Interestingly, the even-odd effect 



.(even) 

^00,0 - - m 



= -ln( 
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is enhanced in the W = (0, 0) winding sector, in agreement with CFT. We also observe that 
the odd curves in both sectors almost coincide on the lattice, as well as in the scaling limit. 
This coincidence can be explained from the CFT result Sit y = £y/ Ly = 1/2: 

st''\y = 1/2) - st'SKy = 1/2) = 2q' + 0{q*) , q = e-\ (37) 
This is a consequence of cancellations at the free fermion point a = 2; in general 

st'^\y = l/2)-stf{y = 1/2) = 2(gV"-g«/4)-2(g2/«-g"/2) + (8/3)(g3/"-g3"/4) + 0(g«) . 
5. Correlations in the SU{N) RVB wave function 

We now start to explore the question of universality, with the motivation of understanding 
whether the results obtained for the dimer RVB state can be adapted to describe the SU (N) 
RVB states. In this section we use quantum Monte Carlo techniques to find the spin-spin 
and dimer-dimer correlators for = 2, 3, 4, 5, generalizing some of the results of [5, 6]. This 
shows that there is no evidence for a transition as is varied, a strong piece of evidence for 
universality. It also allows us to make contact with the results of [21], where an expansion 
around the dimer RVB state is developed. 

5.1. Monte Carlo methods 

To study SU{N)-I{VB wavefunctions, we use an unbiased Monte Carlo technique that 
computes expectation values in an equal-amplitude superposition of states using a Metropolis 
importance sampling scheme [6]. Numerically, the RVB wavefunction is represented in a 
combined VB-spin (or dimer-spin) basis, \Vc) = \Dc)\Zc)^ where \Dc) is a list of site-pairs, 
[z, j], specifying the ends of the valence bonds, and \Zc) is the spin state for each site on the 
lattice. 

In order to obtain expectation values (as described below), two RVB wavefunctions 
are sampled simultaneously, {Vc\ and \Vcf). The two VB states {Dc\ and \D[^) can initially 
be chosen at random; spin states compatible with a non-zero overlap (VclV^/) are then 
constructed. Since the spin basis is orthogonal, we set \Zc) = {Zl^) to ensure a non-zero 
overlap. Then, the Monte Carlo sampling scheme proceeds in two distinct steps. 

After initializing the wavefunction in one configuration \Dc)\Zc)^ the first step is to 
sample a new spin state \Zc) = \Zq). This is done at random, subject to the condition that 
the overlap (VclV^/) remains non-zero. The second step is to sample the VB states \Dc) and 
\D[^). This must be done subject to the constraint imposed by the spin states. As in the 
SU{2) case [6, 7], one has the choice of "local" or "non-local" updating schemes to sample all 
possible VB states [56]. In a local update, the simulation chooses two parallel valence bonds 
on a square plaquette and moves each to the empty bond, subject to the condition that the 
spin-basis is compatible with this new configuration: Fig. 6(a). However, if the spin-basis is 
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(a) 




Figure 6. Possible local updates of the valence-bond (VB) configuration \Dc) in a Monte 
Carlo simulation, where thick lines are the VB states. Figures a) and b) show cases for 
SU{2)^ where in a) the underlying spin configuration (m^ = i^) allows the update shown, 
while in b) it is precluded. In c), an SU{3) case is shown with an additional spin state 
rrii = (top two sites), which precludes the update shown. 



not compatible with this proposed change, as in Fig. 6(b-c) where the right pointing arrow 
represents a different flavor spin for an SU{N > 2) system (e.g. = for SU{3))^ the 
move is rejected, and a new "plaquette flip" is proposed by the simulation. While simple 
to perform, this move is not ergodic in that it cannot change the winding number of the 
system. Therefore, this update is most useful if one wants to fix the winding number to 
do measurements (e.g. for the correlation measurements below, we use the winding number 
W = (0, 0)). In contrast, a non-local loop update moves through the ensemble of states by 
creating a defect at some lattice site and propagating it through the system (subject to the 
spin constraints) until the defect reaches its initial point, thus closing the loop [56]. 

The important point to note is that the Monte Carlo update for the RVB wavefunction 
proceeds in two steps: sampling the spin state \Zc)^ and sampling the VB states \Dc)- Each 
basis configuration is constrained by the other such that the valence-bond remains a quantum 
SU{N) singlet. In contrast, for the classical dimer model, no spin state exists, and dimer 
wavefunctions are orthogonal. 

We now discuss the methods we have used to generalize the SU{2) updates used in 
Ref. [6, 7] to account for the complexities introduced by the SU{N) RVB wavefunctions. 
We do this by considering the similarities and the differences between the SU (2) and SU (3) 
wavefunctions. Given Eq. (9), the singlets on a pair of sites [i^j] can be written as 

Mn=2 = ^[\t)U)j - i;).it)i], t,i= ±1 (38) 
M^=3 =^[\mi)j + \ih\t)j-mo)j], t,;=±i. (39) 

While sampling the spins, there are N different spin flavors to choose for a site on sublattice 
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(a) SU(2) 




Figure 7. Simple transition graphs for RVB wavefunctions. In a), an SU (2) state illustrates 
two different types of transition graph, with "loops" labeled 1 and 2, based on whether a 
spin is +^ or — ^ on sublattice A. VB updates must occur between loops of the same 
type (as in Fig. 6a)). In b), an SU{?>) state has three loop types, depending on whether 
is +1, -1 or on sublattice A. 



A. Once this spin is chosen, all other spins in the same loop are determined by Eq. 39. For 
instance, consider the SU{2) wavefunction in Fig. 7(a). There are two types of loops: one 
with t (loop 1) and another with | (loop 2) on sublattice A. Once the spin basis is chosen in 
the transition graphs, the VB configurations in the bra and ket are independently updated 
as prescribed above according to the same spin basis. Updates, however, are only possible 
between loops of the same type (1 or 2). Therefore, instead of keeping track of the spins on 
each site, it suffices to keep track of the loop labels at each site. 

Then, generalizing this a little further, consider Fig. 7(b). One can easily see that there 
are three loop types: ^^id on sublattice A. So, for an SU{N) system, there are 

different loop types, and we quickly run into ergodicity issues as updates are more readily 
rejected. We find that as A > 4, the number of updates that one can successfully make 
becomes much smaller, making the statistical error bars much slower to converge in Monte 
Carlo. As shown by our data, some estimators for this large N show signs of ergodicity loss; 
this might be cured by a more sophisticated sampling scheme. 

In our VB basis Monte Carlo simulations, we begin by computing the spin-spin and the 
dimer-dimer correlations. While the two-point spin correlation is well-known [57, 58], the 
four-point spin interaction was recently derived in [59]. The expectation values were then 
generalized in [60] for SU{N) spins. Using the notation {ij)L for two sites i^j in the same 
loop in the transition graph and {i)L{j)L for sites not in the same loop [6], 

{Vc\S..S,\Va) _ f i^UJU .^q. 

{Vc\Va) ~\e,,S{S+l) (tj), ■ ^^^^ 

Here, eij = 1 if are on the same sublattice and eij = —1 if they are on different 
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sublattices. We also measure parallel nearest neighbor dimer-dimer correlations, where four 
sites i, fc, / are related hy k = i + I = j + Ga and j = i + mci^^ where a, /3 = x, ^, a 7^ /3, 
and mG{0,l,...}. This simplifies the estimator, which is as follows 

r (r r)- (^c|(S.-S.+eJ(S,-S,+eJ|l^C^) .... 

^dd{ri, rj) - ^^-p^^-^ (41) 

1^2(5 + 1)2 ^ (^^•)^((^ + e„)(i + e„))i or + e„))i((i + e„)i)i 

and otherwise. 

In the next section we present data for these correlation functions for SU{N) with 
= 2, 3, 4, 5, and fit the putative universal exponents associated with their decay. 



5.2. Equal-time correlators 

The connected dimer-dimer correlation function has been recently shown to decay 
algebraically in the SU{2) case in Refs. [5] and [6]: 

Qd(ri,r2)~^ ^ (43) 

1^1 — V2\ 

at distances large compared to the lattice spacing, but small compared to the torus size 
L = Lx = Ly. Using several diflFerent methods, the exponent was found to be approximately 
a 1.2, independent of the direction. We extend some of these results to the SU{N) case, 
mainly focusing on the correlation between two parallel dimers in the transverse direction 
(ri — r2 = ix^x)- Out quantum Monte Carlo results are shown in Fig. 8 for a L x L torus, 
with L = 64. The data clearly indicate that the algebraic decay remains for = 2, 3, 4, 5, 
with a increasing with A^. Since SU (N) dimer correlations reduce to pure dimer correlations 
in the limit ^ oc, we know that limA^^cx) o^at = o^d = 2, and our data is consistent with the 
smooth interpolation between N = 2 and oc. For each A^ we extract the leading exponent 
by a fit (solid blue lines in Fig. 8) to Eq. (43) in the regime 1 <C <^ and for odd ix- We 
find a2 = 1.21(7), as = 1.43(7), = 1.54(8), = 1.67(9). Notice that as A^ gets larger, it 
becomes more and more difficult to converge the Monte-Carlo data in a finite time, resulting 
in a loss of precision on the corresponding exponent. 

We also observe the presence of a dipolar term which behaves as (—1)^^^"^ for all A^. This 
direction-dependent contribution is well known in the dimer model [24] , and also appears for 
the SU{2) RVB state [6]. In the dimer model and for the transverse correlations considered 
here, it cancels exactly the universal term (43) at even distances, so that the dimer-dimer 
correlation behaves as [10]. In the SU{N) RVB case the exponent a is different and this 
does not happen. The dipolar is then a subleading term, which introduces finite-size effect 
when trying to extract the exponents. For example we get a2 — 1.14 at even distances, 
which is slightly smaller than the previous computed value 2::^ 1.20 obtained for larger system 
sizes [5, 6]. 
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Figure 8. Dimer-Dimer correlations on the 64 x 64 torus for N = 2,3,4,5. Dashed blue 
lines is a linear fit to Eq. (43), in the "linear regime" ix — 5, . . . , 15, while the red solid lines 
are a fit to the conformal scaling (Eq. (44)) in the same region. 



These values are in good agreement with the expansion around the dimer RVB state 
developed in [21]. This is a cluster expansion of the loop model, describing the non-trivial 
inner product in (11) perturbatively in 1/A^ in terms of a classical interacting dimer model. 
To first nontrivial order, this model is exactly the one studied in [23, 24], and so the Monte 
Carlo results in Fig. 26 of [24] can be utilized; in their notation X2 = and W = 1 + 1/N. 
This gives 0^3 1.4, 1.52, and 1.6, in addition to the already reported 0^2 — 1-22. 
Our numerical results therefore confirm that the leading-order term in this expansion gives 
accurate results. Note that higher order terms could slightly increase these values, as is also 
shown in the supplementary material of [21]. 

A presumably more- accurate value for a can be found relaxing the constraint i ^ 
and fitting the data to the full curve predicted by conformal field theory, not just the straight- 
line segment. While the two-point function on a torus is not fixed uniquely by conformal 
invariance, the appropriate form is known for the free-boson field theory for any [47]. Thus 
if we assume that the free-boson/Coulomb gas results still apply with a different = 1/a, 
the SU (2) dimer-dimer correlators are modified to 

I^||(4)-/(4/I^x,r)-", (44) 

with ix ^ 1. / is a universal function of the two dimensionless ratios ix/Lx and r = iLy/Lx 
and can be expressed in terms of a Jacobi theta function defined in Appendix B. We have 

fiijL,, t) = sin [{2n + l)nijL,] (45) 

n=0 
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Notice that while this function reduces to f{£x/Lx^r) = 8m{7T£x/ Lx) in the hmit of a thin 
torus {Ly/Lx ^ 1), the approximation remains exceUent for Ly/Lx = 1. Performing the 
fits to Eq. (44) in the same region 1 <C <^ as before (red sohd curve in Fig. 8), we 
observe that the data reproduces very weU the free-boson CFT result, even the upturn when 
£x is of order Lx- The fact that the CFT result seems to still hold here provides additional 
evidence in favor of universality. However, the exponents determined this way are slightly 
diflFerent; for example we find 0^2 — 1-26 in the odd case, and a2 — 1.18 in the even case. 
This small discrepancy with the previous SU{2) results [5, 6], also suggested by the large 
calculations [21], could possibly be resolved by studying larger systems. 
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Figure 9. Spin- spin correlation function on a 64 x 64 torus. We look at spins separated by 
a distance along the x— direction: ri — r2 = tx^x- 

We also studied spin correlations, which as discussed in section 2 are known to decay 
exponentially in the SU{2) case: {S{0)S{£x)) ^ exp(— 4/<^) [28, 5]. We find the same 
behavior in the SU{N) case, with our results displayed in Fig. 9. We also observe that ^ 
decreases with A^, compatible with the intuition that the RVB state becomes more and 
more spin-disordered as increases, approaching zero correlation length in the dimer 
limit ^ oc. Fitting the data gives ^{N = 2) = 1.30(5) (compatible with [5]) and 
^(A^ = 3) = 0.68(5), ^{N = 4) = 0.56(3) and ^{N = 5) = 0.46(2). This lack of long-range 
Neel order is consistent with quantum spin-liquid behavior. 

6. Entanglement in the SU{N) RVB wave function 

Throughout this paper we have emphasized the similarities between the SU (N) and dimer 
RVB states. These make it plausible that the two-cylinder Renyi entropy for the SU (2) RVB 
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state in the scaling limit can be obtained from the free-boson result (16). If this is true, 
this is a much stronger statement than the results for the critical exponent a as a function 
of in the previous section. The reason is that even though correlators in the classical 
dimer model can be obtained by taking the ^ oc limit of the loop model defined by the 
SU{N) RVB inner product, the Hilbert space of the two is completely different; the latter 
being much larger. Thus a dimer crossing the boundary would carry much more entropy. 
Moreover, the reduced density matrix for SU (N) has blocks with non-zero total spin, while 
these are by construction forbidden in the quantum dimer model. Thus it is very unlikely 
that non-universal quantities will be the same in the two cases. 

All our results are consistent with the subleading part of the two-cylinder Renyi entropy 
being universal. Thus it is possible that in the scaling limit, this piece still could be given 
by the conformal field theory result, with the identification of k coming from the correlation 
length results. In this section we provide evidence for such universality by finding the curves 
for arbitrary a, and extending previous numerical results [7] for SU{2) to the SU{N) case. 

6.1. Monte Carlo evaluation of the Renyi entropies 

In order to calculate the Renyi entropies, the Monte Carlo sampling algorithm of Section 5.1 
must be modified to incorporate a replica incA; procedure, first described in Ref. [35]. Namely, 
in order to calculate the n-th order Renyi entropy Eq. 1, one measures the expectation value 
of a Swap A (or permutation) operator acting between n copies of the system. For the SU (N) 
wavefunctions studied here, we restrict ourselves to n = 2. 

A naive sampling of the Swap a operator is straightforward. One simply identifies the 
spatial region A and its complement 5, then applies the operator to a sampled RVB state, 
SwapA\Vc) = \Vc*) (where each state \Vc) contains n = 2 replicas). The second Renyi 
entropy can be calculated from [35] 

p\ = ^-^^^ = iV"=--"' , (46) 

where ni is the number of loops in the transition graph of the two un-swapped RVB states 
(the denominator) , ngwap is the number of loops in the transition graph using one swapped 
state, and A^ denotes the number of spin fiavors from SU{N). Thus, the reconfiguration 
of the valence-bond endpoints is the important ingredient of the SwapA operator for this 
measurement. Note, since the Swap a operator serves to swap all basis states in the spatial 
region A, spin states \Zc) can eflFectively be ignored when performing this measurement; they 
are only used as a condition to sample the original un-swapped VB configurations. 

However, one discovers quickly as the lattice size (and particularly the size of region ^4) 
grows, sampling statistics become exponentially poor when using the naive Swap a operator 
to calculate 5^2. One must therefore employ some variant of the ratio trick described in 
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Ref. [35]. As discussed in detail there, this involves sampling 
{Vc\SwapA\Vn) 

{Vc\Swapx\Vny ^ ^ 

where X is a spatial region smaller than i.e. X C A. In order to implement this 
measurement in the RVB wavefunction, a completely different simulation must be run with 
the overlap (yc\Swapx\Vn) occurring in the weight. 

The essential point of the ratio trick then is that simulations must be run with the 
basis states corresponding to region X already swapped before measurement of the Swap a 
operator. Since we work in a combined VB-spin basis, one imagines the simulation beginning 
with two states, \Vc) = \Dc)\Zc) and \Vn') = |i^7^/)|.^7^/). However, the simulation is 
performed with (yc\Swapx\Vn) in the weight, meaning that \Dn')\Zn') = Swapx\Dn)\Zn) - 
Therefore, it is important to sample spin states such that the overlap {Zc\Zn') is nonzero. 

Note that, when performing the Monte Carlo update on the valence-bond state |i^7^/), 
one does not have to consider the overlap of (yc\Vn') since the overcompleteness condition 
ensures that the inner product of new valence-bond configurations (with the same spin 
configuration) are always non-zero. Hence, it may be convenient to consider the "un- 
swapped" state \Vn) for the purposes of sampling the VB configuration, where one can 
straightforwardly perform local (plaquette-fiip) or non-local (loop) updates on the VB state 
as described in Section 5.1. 

We employ ratio-trick evaluations of the Swap a operator to calculate the scaling of the 
second Renyi entropy in the SU (N) RVB wavefunction. In the next subsection, we compare 
our results from these Monte Carlo simulations to those obtained from conformal field theory. 

6.2. Comparison of numerical results and CFT 

Here we compare our quantum Monte Carlo results with curves found by using Eq. (16) 
with partition functions coming from conformal field theory for arbitrary stiffness = 1/a. 
This two-cylinder entropy is that of the quantum Lifshitz state, the continuum limit of 
the quantum six-vertex state. If the universality discussed above holds true, they will also 
apply to the continuum limit of the square-lattice SU{N) RVB states, with appropriate 
identification of a in terms of A^. 

Let us begin with the SU (2) case, and present our results for 5^2 both for unrestricted 
windings and for the W = (0, 0) winding sector. The quantum Monte Carlo data is shown 
in Fig. 10. We observe a strong similarity to the dimer data. In particular, there is a strong 
even-odd effect that becomes even bigger in the W = (0, 0) sector. This suggests 5'2 for the 
SU{2) RVB state lies in a locked phase, just as its dimer counterpart does. This scenario 
is also supported by exact diagonalizations on small systems, which show that the biggest 
eigenvalue of the reduced density matrix is non-degenerate: /S.su{2) = — ln(]9i/]9max) ^ 2.67 
for a 4 X 4 torus cut in two halves, and therefore dmin = 1- (The exact value for dimers is 
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Figure 10. Numerical extraction of the universal shape of 5*2 as a function oiy = ^y/ Ly for 
the SU{2) RVB wave function and L = = Ly = 12, 18, 24. Red and blue curves are the 
CFT results, using Eqs. (48,49,50,51). Left: W = all sector. Right: W = (0,0) sector. We 
observe in this case a stronger even-odd effect on the lattice, consistent with the expected 
result in the scaling limit. 



Adimers = 2 In TT 2.29 in the scaling limit [61].) Such a large entanglement gap has also been 
observed recently for the SU (2) RVB state, in the closely related infinite-cylinder geometry 
[62] . We also refer to [62] for a detailed study of the entanglement spectrum in this geometry. 

With no restrictions on the winding numbers, the conformal field theory result for the 
two-cylinder Renyi entropy at general a is 
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We have made the assumption that the boundary conditions are the same as those for dimers 
(i.e. a = 0, 1/2 for even and odd respectively for all a). 

The curves (48) and (49) with a = 1.2 are plotted in fig. 10 along with the SU{2) Monte 
Carlo data. The trend in the data is clearly toward these curves, but as with the dimer case, 
the finite-size effects are very strong. Several additional interesting features become apparent 
by plotting the CFT curves (48) and (49) at a = 2 and a = 1.2, corresponding to the dimer 
and SU{2) RVB exponents respectively. As is apparent from fig. 11 (left), the curve for the 
even branch is essentially independent of the a over a large range, and resembles very much 
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Figure 11. CFT predictions for the universal shape of S'2, as a function oi y = iy/Ly.Left: 
Red sohd hues correspond to the dimer RVB state with a = 2, while the blue dashed lines 
are the SU (2) RVB case, with a = 1.2. We also plot in black the usual Id result a-\-b In sin Try 
with appropriate a and b. Right: Prediction for a = 1.2. Blue dashed lines is the ly = (0, 0) 
sector, while red solid lines show the prediction with no restriction on winding numbers. 

the Id result 5^2 oc lnsin(7r^). The shape of the odd curve only weakly depends on a, but the 
curve does appreciably shift in this range. We also observe -fig. 11 (right)- that contrary to 
the dimer case (a = 2), the two odd curves corresponding to a = 1.2 are not nearly identical 
in the W = (0, 0) and W = all sectors. This feature is also in good qualitative agreement 
with the data shown in fig. 10. 

Also plotted in fig. 11 is the curve a + 61n(sin7r^), with a and b fitted to make the curves 
as close as possible. This curve is the Id Renyi entropy for a conformal field theory with space 
a circle of circumference L split into segments of length yL and (1 — y)L. The logarithmic 
term is the leading term in one dimension, with its coefficient b proportional to the central 
charge [55]. One might expect that our 2d curve would resembles the Id curve when the 2d 
system is obtained by coupling together critical Id systems maintaining the criticality. This 
for example is the case for fermions with tt flux discussed in [7]. However, the square-lattice 
quantum dimer model is gapless only in the two-dimensional limit (just like the Goldstone 
modes for the Heisenberg model also discussed in [7]), so we see no particular reason for the 
curves to be so similar. Nonetheless, as is obvious from the figure, the two curves are quite 
close, although clearly not identical. 

The same picture still seems to hold for the SU{N) RVB state. We plot S2 data for 
L = 16 and S'[/(2,3,4) is shown in Fig. 12. As can be seen, the error bars get bigger and 
bigger for larger A^, but we still observe a clear even-odd effect. This strongly suggests that 
we are still in a locked phase, and that Eq. (48) applies in the thermodynamic limit. 

Extracting an accurate exponent from the 5^2 data is difficult, due to the combined 
finite-size effects and statistical errors. A possible try would be to look at the even-odd 



29 



10 


■ 




l'= 16, S'f7(2) 

L = 16, S'[/(3) HHe-H 

L = 16, S'[/(4) H-^H 




■ 


9 

8 


\ i, 

■. \ / (f \ , 


A 


A A 

V T V 


-A 


, / 05 \ / / 


7 


- \\ // \V 


9. 


9, 




V / w // - 


6 


• i : 




V * V 




* A ^- 


5 


\ 6 / \ 








/ \ d / 


4 













2 4 6 8 10 12 14 16 



Figure 12. Numerical data for the second Renyi entropy 5*2 in the SU (2, 3, 4) case. As can 
be seen, even on relatively small system sizes the numerical data is not totally converged. 

difference for y = 1/2. On the lattice it is defined as Sn{Lx^ Ly) = \Sn{Ly/2 + 1) — Sn{Ly/2)\ 
for even Ly. This diflFerence should go to a positive constant in the thermodynamic limit. 
From Eq. (48) we predict 

.„(r,„)^^ infill (52) 



l-n \e^{aT/2) 

Using our biggest system-size even-odd data in (52) yields a 2.27 for dimers and a 1.57 
for SU{2) RVB. These numbers can be improved by extrapolation, but the strong finite-size 
effects described above prevent an accurate estimate of these numbers. 

With such strong finite-size effects, we cannot make a definitive conclusion. Our results, 
however, provide qualitative evidence for Eq. (16) applying to the SU{2) RVB state. They 
also support the convention of [7] that the subleading piece of the two-cylinder Renyi entropy 
is universal. 

Unfortunately, we do not know how to derive these results as we did for the dimers. 
The mapping of Ref. [15] to a classical Renyi-Shannon entropy is not valid for the SU{N) 
RVB states; because of the non-orthogonality of different valence-bond configurations, the 
Schmidt decomposition of the RVB state is much more complicated. Moreover, a naive 
extension of the result of [18] would then yield a critical value ric = a/2 0.6 for SU{2). 
This would imply an even-odd effect even for the von Neumann entropy, in violation of strong 
subadditivity [19]. The formula Uc = <^min/(2^) = <^min^/2 therefore cannot hold in general 
for models that cannot be mapped onto a classical Renyi-Shannon entropy. The fact that it 
does apply to quantum six-vertex states does not contradict strong subadditivity, because 
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in this case rfmin = 2 and a > 1, so nc > 1 in this case. § 

A possible way to test our CFT results while sidestepping this difficulty would be to 
look at another quantity, defined as follows. Consider another SU{2) RVB wave function 
1^') living on the torus, but where all valence bond states \Vc) with singlets crossing the 
boundaries between A and B have been removed. Such a wave function can be rewritten as 
the tensor product of two RVB- wave functions living respectively on cylinders A and B: 

\^') = \^{A))0\^{B)) . (53) 

Consider then the logarithmic (bipartite) fidelity [63] defined by taking the scalar product 
with the original (torus) RVB wave function: 

F = -ln(|(^|^Ol')- (54) 

This fidelity enjoys several properties very similar to the entanglement entropy, like a generic 
area law, and a universal logarithmic divergence for one-dimensional quantum critical points. 
Using Eq. (8), it can be expressed using classical partition functions for the RVB loop gas. 
In this case loops longer than the dimers can cross the boundary between A and 5, but 
since the dimers dominate, we expect our main result (16) to apply to this quantity, upon 
removing the n/{n — 1) factor (or equivalently setting n = oc). Using the arguments of 
section 3.2, one can also check that the equality 

F = 5oo (55) 

holds exactly on the lattice for the dimer RVB state. While this is not true anymore for 
SU{N) RVB states, we expect (55) to remain true at the level of universality. Therefore, 
Eqs. (48,49) should describe the shape of both Soo and F, when varying the respective sizes 
of the cylinders A and B. Studying this fidelity, for example using quantum Monte Carlo 
methods similar to the ones presented here, could probably shed some more light on the 
phase transition scenario. 

7. Conclusion 

In this paper we studied a family of resonating-valence-bond (RVB) states on the square 
lattice, focusing in particular on the Renyi entropy Sn for dividing a torus into two cylinders. 
This provides an elegant probe into the physics of such gapless states. We found several 
interesting characteristics of their quantum entanglement, by deriving exact results for the 
dimer and quantum Lifshitz RVB states, and by numerical work on the SU{2) and the 
SU{N) RVB states on the square lattice. 

One of our main results is confirmation that the Renyi entropy for RVB states exhibits a 
two-branch structure, when the length of the entangled region varies between an even and odd 

§ Notice that naively applying ric = to an RK wave function built out of interacting dimers [23, 24] 

would, as with the SU{N) RVB, violate strong subadditivity. 
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number of lattice sites. We computed this exactly for the dimer RVB functions by mapping 
the Renyi entanglement entropy to a classically-computable Renyi-Shannon entropy, and 
demonstrated that the two-branch structure exists for n > 1 (but not for n < 1). By 
measuring the Monte Carlo expectation value of a "swap" operator on a replicated lattice, 
we find 5^2 for the SU {N) RVB state for = 2, 3, 4, and find the same two-branch structure. 
It would be interesting to extend the Monte Carlo calculations to higher n, where the replica 
trick can still be used, and to study the fidelity of Ref. [63] in this context. 

The two-branch shape of the Renyi entropies has an interesting consequence: strong 
subadditivity is violated for the "odd" branch, which is a concave-up function of the linear 
size of the entangled region. We demonstrate that this is far from a numerical artefact even 
away from the dimer case, by calculating the shape of these branches in an effective conformal 
field theory. This effective theory involves partition functions that can be evaluated exactly 
in the continuum, and exhibit the same two-branch structure observed in lattice calculations. 
This emphasizes the universality of the result and confirms that the two-branch Renyi 
entropies are not the result of finite-size lattice effects, contrary to what generically happens 
in one-dimensional systems [64, 65]. 

From this analysis, we have elucidated two important points regarding the shape of the 
Renyi entropy as a function of linear size of the entangled region. First, the two-branch 
structure does not occur for every Renyi index n, but only for values strictly larger than 
some (non-universal) nc, which must be greater than unity since the von Neumann entropy 
Si is strictly constrained by strong-subadditivity (while the n> 1 Renyi entropies are not). 
We believe that in the square-lattice RVB, 1 < nc < 2; since our Monte Carlo methods 
require n > 2, they will always find a two-branch structure as observed here. It would be 
useful to change rfmin by studying the RVB wave function on other bipartite lattices. For 
example on the honeycomb lattice rfmin = 3 (and therefore nc = 3^ = 9) for dimers, and 5^2 
for RVB would likely lie in the replica phase, and would not give rise to an even-odd effect. 
This prediction should be tested by future Monte Carlo studies. 

Second, we have found an exact expression for the shape of the Renyi entropy curves 
in a geometry where a linear size of the entangled regions grows as I. The fact that it is 
consistent with our numerical data for the SU (2) and dimer RVB wavefunctions is a strong 
indication that it is universal. We see no particular reason why this universality should be 
a characteristic only of RVB wave functions, so this conjecture should be tested on future 
gapless wavefunctions. Indeed, DMRG results for the 2d transverse-field Ising model critical 
point have already exhibited the same type of behavior [66] . As observed in previous studies 
of gapless systems [7], these curves are well approximated by the celebrated universal one- 
dimensional result ln(sin(7r^/L)); however, deviations were observed. Our exact 2d results 
from effective field theory indeed confirm that the Id result is only approximate in the RVB 
case, but also that it is a very good approximation. 

We have also explored the crossover between the SU{2) and dimer RVB states by 
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studying the SU (N) states. For equal-time correlators, it is possible to interpolate between 
the SU{2) and dimer states by taking ^ oo, but there is no guarantee that there is no 
transition at some value of A^. Previous studies demonstrated that the SU{2) RVB state 
exhibited algebraic decay with an exponent a2 ^ 1.2, while the dimer RVB has a^o = 2. We 
showed in this work that this exponent evolves smoothly as is varied, with no evidence for a 
transition. This is strong evidence of universality, allowing us to link various physical aspects 
of the SU{2) case with exact results on the dimer RVB. Perhaps even more importantly, we 
also showed that this correspondence at least qualitatively holds for the putatively universal 
part of the two-cylinder Renyi entropy as well. This to us is strong evidence both for the 
universality and the effectiveness of using dimer models to probe at least some highly non- 
trivial aspects of the SU (2) RVB state. 

It is worth recalling that the SU{N) and dimer RVB states are the simplest 
wavefunctions exhibiting features of gapless spin-liquid behavior. Since gapless spin liquids 
are long-range entangled states not amenable to classification by simple quantities such as 
the topological entanglement entropy, we hope that some of the universal scaling properties 
uncovered here will be of use in future efforts to classify and characterize them in theory, 
numerics, and experiment. 
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Appendix A. Renyi- Shannon entropy for dimers on the torus 



We wish to compute the foUowing classical entropy 



Sn — 
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n 



(A.l) 



(7,11 



for any real n. The probabilities p^r,^ we are interested in are given by a ratio of partition 
functions, which can handily be expressed using a transfer matrix 

Z Tr T^y ^^'^^ 

T is the transfer matrix of the dimer model, and acts on the vector space generated by dimer 
occupancies on vertical edges along a horizontal line: {a\T\b) = 1 if configuration \a) and \h) 
are compatible, otherwise. The denominator is calculated in Appendix C. Each factor on 
the numerator of Eq. (A. 2) can be evaluated using free fermions techniques [67, 24, 15]. 

The mapping goes as follows. Let us first consider a particular configuration 
with staggered dimers, which we call reference configuration (see fig. Al on the left). 
Superimposing the reference on any other dimer configuration generates a collection of non- 
intersecting lattice paths, represented by black lines in fig. Al (right). 
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Figure Al. Illustration of the mapping onto free fermions for open boundary conditions. 
Left: reference configuration with staggered dimers. Right: transition graph generated by 
superimposing the reference configuration on a given dimer configuration (in thick blue). 
This generates a collection of 4 non-intersecting lattice paths going upwards, which may be 
expressed in terms of fermions (black zigzag lines). Here |<j) = C2c\cqcI\0). 

There is a one to one correspondence between the dimer configurations and the lattice 
paths, provided the latter obey a certain set of rules. To clarify them, it is most convenient 
to reformulate everything in terms of fermions, because their statistics will naturally encode 
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the non-intersection constraint. Define a fermion as a vertical link occupied by a reference 
(red) dimer only, or a real (blue) dimer only. These are represented by black zigzag lines in 
fig. Al. I a) and |//) can be rewritten using fermion creation operators as 

\a) = |xi,X2,...,X7v) =4i42---4^|0) (A.3) 

I//) = |2/i,2/2,...,?/7v) = 4jA2"'4jn\^) (^-4) 

where the Xi and yj label the ordered positions of the fermions. Introducing of shift of +1 
lattice spacing from one row to the other, it is easy to check from fig. Al that the transfer 
matrix satisfies 

T|0) = |0) (A.5) 
Tcl^_,T-' = (A.6) 

^4j^ ^ = 4j + 4j+i + 4j+2 i^-^) 

Eqs. (A. 5, A. 6, A. 7) conserve the number of fermions, and encode the dimer close-packed and 
hard-core constraints on the square lattice. To account for periodic boundary conditions 
we also have to impose 4^,+! ~ {—l)^^^c[> Rewriting this in matrix form using Einstein's 
summation convention 

Tc\T-^ = M,jc], (A.8) 

we successively get 

T^T-'={M')^^c] (A.9) 

and 

{a\T%) = {0\c,,c,, . . . c^^T'clcl . . . 4|0) (A.IO) 

= (0|c.,c., . ..c,,{M%,,,cliM%,,,cl . . . {M%,,,cl\0) (A.ll) 

By applying the Wick's theorem and using (0|ca;C^|0) = 5xy^ we finally get 

{a\T%) = Jet^ [{M%^,^] (A.12) 

This type of result is referred to as the Karlin-McGregor[68] or Lindstrom-Gessel-Viennot[69, 
70] lemma in the mathematical literature. In the end, Eq. (A. 2) reduces to a product of two 
determinants, which may be computed numerically. 
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Appendix B. Dedekind-eta and Jacobi-Theta functions 

For a given modulus r, we introduce the squared nome: 

q = e^^^^ (B.l) 
The Dedekind eta function is defined as 

oo 

k=l 

and the four Jacobi Theta functions are given by 

e,{z\T) = J^(_l)-V2^i(n+l/2)^g(2n+l)i. ^33) 

nGZ 

00 

= 2r){T)q^/^ sin ^ JJ [1 - 2 cos(2^)g" + g^"] (B.4) 

n=l 

e^{z\T) = ^g|(-+V2)^e(2n+l)^. (35) 

00 

= 2r]{T)q^'^ cos z JJ [1 + 2 cos(2z)g^ + q^""] (B.6) 

n=l 

e3(z|r) = ^qh\^^^^ (B.7) 

nGZ 

00 

= r]{T)q-^'^^ n [1 + 2 cos(2z)g"-^/' + g'"-^] (B.8) 

n=l 

^M^) = ^{-iTq'^'^'e^'''' (B.9) 

00 
n=l 

To simplify the notations a bit, we also set 

e^ir) = e,{0\T) , u = 1,2,3,4 (B.ll) 

These functions obey the following nice "modular" transformation properties 

^i(r) = -i(-«T)-^/2^i(-l/r) (B.12) 

0^^r) = {-tT)-'/'e,{-l/T) (B.13) 

0,(r) = {-tT)-'/'e,{-l/T) (B.14) 

^4(r) = (-ir)-V2^2(-l/r) (B.15) 

ri{T) = (-ir)-i/2^(-l/r) (B.16) 
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Appendix C. Derivation of CFT partition function from the lattice dimer model 



Appendix C.l. Torus 

The goal of this section is to recover the CFT partition function on the torus, using the exact 
solution for dimers in terms of transfer matrix, and keeping track of the winding numbers 
Wx and Wy. To do that we need to solve a slightly more complicated problem than that 
in Appendix A, introducing a fugacity h (resp. h~^) for horizontal dimers whose left site 
belongs to the even (resp. odd) sublattice. We will also suppose and Ly to be even. With 
this at hand the transfer matrix now satisfies 

m = io) (C.l) 

Tcl^_,T-' = 4^. (C.2) 
Tcl^T-' = b 4 + 4^.+i + b-'cl^^, (C.3) 

The 6's provide information about Wx^ while the number of fermions does that for Wy. 
Defining as before T(^-T~^ = MijCj^ diagonalizing T amounts to diagonalizing M, which can 
be done in Fourier space, carefully taking into account the boundary condition on fermions 

cL+i = (-l)^+M , N = J^c]cj (C.4) 

i=i 

We write the number of dimer coverings in the {Wx,Wy) sector as Zw^,Wy- Setting 
b = e-'^"^/-^'^^ the winding generating function is given by 

ZtoU^..Uy)= 5] ZH^,,H^^e-^(^-^-+^^^^) (C.5) 

= Tr |^e~^^^^^~^^/^^/^T^^ (C6) 

If we now denote by the set of fermionic operators that diagonalize the one-particle 
transfer matrix 

TdlT-' = A,4, (C.7) 
the Xk are then the eigenvalues of M. Using this, the transfer matrix can be expressed as 

T = l[{l + [Xk- 1]44) (C.8) 
k 

Plugging this expression in Eq. (C.6) we get|| 

ZtoUu., Uy) = \ {Z+ + Zf ) + ^ {Z+ - Zo") (C.9) 

II One can check easily that expanding simultaneously the first two terms in Eq. (C.9) generates all the 
eigenvalues of the transfer matrix in the even fermion sector. The same goes for the last two terms and the 
odd-fermions sector. 
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with 



and 



Afc = cos k + \/l + cos^ 
(2m + + iMy)7r 



m 



-LJ2,...,LJ2-1 



(C.IO) 

(C.ll) 
(C.12) 



It is also possible to rewrite as 



m=0 



e^^- + e-^^- ± 2rr icos 



(2m + V + iUy)Ti 



(C.13) 



where T^(x) = cos(marccosx) is the m— th Chebyshev polynomial of the first kind. This 
expression allows for convenient numerical evaluations on large system sizes. 

We now wish to extract the universal CFT partition function i.e the constant term 
in the asymptotic expansion of Z: 

^torus(^., Uy) - A^-^-S^-C^-Ztorus(^x, Uy) (C.14) 

A possible way to do so would be to combine Euler-Maclaurin expansions and maj oration 
techniques, as has already been done for the honeycomb lattice[71]. Here we obtain 2^ in a 
more heuristic (and non rigorous) manner, noticing that universal properties have to come 
from low energy excitation near the Fermi momenta kp = ±7r/2. Linearizing ln(A/e) around 
the two kp and extending the products in C.9 over all integers, allows to go get after a long 
calculation, very similar to that in [71]: 

e^{i^uj2\f/2)e^{i^Uy/2\T/2) 



{2'^mrf'^ r]{T) 



(C.15) 



where r iLy/Lx and f —1/^ iLx/Ly. 



Appendix C.2. Cylinder 

The cylinder geometry is slightly simpler than its torus counterpart, because only the 
windings along y are allowed. To still express the partition function as a trace, we look at 
the transfer matrix T acting on the configuration of rows of length Ly with open boundary 
conditions. This way, the number of fermions keeps track of the winding number Wy. The 
winding generating function can be expressed, after some algebra, as 

Ly 

Z,y,{uy) = e-(^^/2)% 11(1 + e^^yfi^^) (C.16) 



m=l 



COS 



+ Wl + C0S2 

Ly + lJ V \Ly + l 



(C.17) 
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Here the Fermi momentum is = 7r/2, and the CFT partition function can be accessed 
using the hnearization trick. We however have to distinguish between the even and odd Ly. 



Even case: It is most convenient to rewrite Z as 



Ly/2 



171=1 



Ly/2 



J] 1 + 2cosh(7r«^)^-^^ + ^-2^^ 



m=l 



(C.18) 



Combining the Euler-Maclaurin formula on the first product, and using the hnearization 
procedure around m = 1/^/2, we get 



(C.19) 



(C.20) 



Zcyiiuy) = f[ [l + 2cosh(7rtiy)e"'^^^"^/^^ + e~^^^^~^^^'^ 
p=i 

which gives, using the infinite product representation of the Jacobi Theta function, 

^3(^7^^t,/2|f/2) 

We recover the Eq. (31), upon setting Uy = and performing a modular transformation. 

Odd case: We use the same method, but care must be taken because of the presence of a 
zero-mode for m = {Ly + l)/2 We have 



{Ly-l)/2 



X 



{Ly-l)/2 



II l + 2cosh(^2i^)/i^-+/i^^- 



m=l 



(C.21) 



(C.22) 



(C.23) 



^cyiK) = 2cosh(7rV2) n 

m=l 

which yields after linearization 

oo 

TcLx -J — w- r '^Lx 

ZcY\{uy) = 2e ^^^y cosh(7r2i^/2) I I U + 2 cosh(7r2i^)e ^ + e 

p=i 

In the end we obtain: 

e2{z7rUy/2\f/2) 

and once again recover the CFT result (32) after setting Uy = and modular transformation. 
Appendix C.3. Zero-winding sectors 

With the winding generating functions at hand, the calculation of 5^(y, r) in the W = (0, 0) 
winding sector becomes straightforward. For example we have in the odd case (recall 
r = iLy/Lx and f = — 1/r) 

02(0 

cyl 



^(4^2(1^)' 



(C.24) 



W={0,0) 

where we have selected the constant term in the product of the two theta functions. This 
allows to recover Eq. (50) and Eq. (51) after once again performing a modular transformation. 
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